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LONG-TERM  GOALS 

Our  long-term  technical  goal  is  to  develop  and  test  the  Tangent  Linear  Model  (TLM),  Representer 
Model  (RPM),  and  Adjoint  Model  (ADM)  versions  of  ROMS  (Regional  Ocean  Modeling  System)  and 
TOMS  (Terrain-following  Ocean  Modeling  System)  for  variational  data  assimilation,  ensemble 
forecasting,  and  stability  analysis.  The  primary  focus  is  to  develop  a  general  platform  for  strong  and 
weak  constraint  4D  Variational  data  assimilation  (4DVar);  to  develop  an  ensemble  prediction 
capability  based  on  optimal  perturbations  and  stochastic  optimals;  and  to  develop  stability  analysis 
tools  based  on  eigenmodes  and  singular  vectors  to  explore  the  role  of  environmental  stochastic  forcing 
in  shaping  ocean  circulation.  Our  long-term  scientific  goal  is  to  model  and  predict  the  mesoscale 
circulation  and  the  ecosystem  response  to  physical  forcing  in  the  various  regions  of  the  World  Ocean 
through  state  estimation. 

OBJECTIVES 

The  objectives  and  scientific  goals  of  the  proposed  research  are: 

1.  To  explore  the  factors  ( e.g .  uncertainties  in  initial  conditions  versus  those  in  surface  forcing  and 
boundary  conditions)  that  limit  the  predictability  of  the  circulation  in  regional  ocean  models  in  a 
variety  of  dynamical  regimes; 

2.  To  develop  state-of-the-art  variational  data  assimilation  platforms  (strong  and  weak  constraint 
4DVar)  and  gain  experience  in  regional  ocean  applications; 

3.  To  develop  ensemble  prediction  techniques  for  regional  ocean  models. 

APPROACH 

This  is  a  collaborative  effort  involving  Dr.  Andrew  M.  Moore  at  University  of  California  at  Santa 
Cruz,  Drs.  Arthur  J.  Miller  and  Bruce  D.  Comuelle  at  the  Scripps  Institution  of  Oceanography,  and  Dr. 
Heman  Arango  at  Rutgers  University.  To  address  the  aforementioned  goals  and  objectives,  we  are 
using  a  newly  developed  suite  of  tools  that  utilize  ROMS/TOMS  tangent  linear  and  adjoint  models. 
These  models  and  tools  were  developed  under  the  support  of  previous  ONR  funding. 

To  address  objective  (1),  we  are  using  the  ideas  of  Generalized  Stability  Theory  (GST)  in  order  to 
identify  the  most  unstable  directions  of  state-space  in  which  errors  and  uncertainties  are  likely  to  grow. 
Specifically,  for  initial  and  boundary  condition  errors  we  compute  the  singular  vectors  of  the  TLM 
propagator,  while  for  surface  forcing  we  compute  the  stochastic  optimals.  By  examining  the  details  and 
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dynamics  of  these  structures  we  can  learn  much  about  the  factors  that  limit  the  predictability  of  the 
circulation. 

To  address  objective  (2),  we  are  using  several  4-dimensional  variational  data  assimilation  schemes  that 
have  been  developed  for  ROMS/TOMS.  For  cases  in  which  the  dynamics  are  imposed  as  a  strong 
constraint  ( i.e .  no  model  error  assumed)  we  use  an  incremental  4DVar  approach  (IS4DVAR)  similar  to 
that  used  operationally  at  some  numerical  weather  prediction  centers.  In  the  case  where  errors  are 
admitted  in  the  model  we  use  an  indirect  representer-based  weak  constraint  4DVar  algorithm 
(W4DVAR)  and  a  weak  constraint  Physical  Space  Analysis  System  (W4DPSAS).  W4DVAR  is  based 
on  the  Oregon  State  University  Inverse  Ocean  Model  (IOM)  of  which  ROMS/TOMS  is  also  a 
component  (Di  Lorenzo  et  al,  2006;  Muccino  et  al,  2006).  The  IOM  requires  an  additional  version  of 
the  model  that  computes  a  finite  amplitude  linear  estimate  of  the  total  state  of  the  system  as  opposed  to 
perturbations  about  some  existing  solution  of  the  nonlinear  ROMS.  This  second  linearized  form  of 
ROMS/TOMS  (denoted  here  as  RPM)  has  been  developed  under  a  separate  NSF  funded  effort. 

To  address  objective  (3),  we  will  use  the  optimal  structures  identified  in  (1)  using  GST  to  construct 
ensembles  of  model  forecasts  following  the  approach  used  operationally  at  some  numerical  weather 
prediction  centers. 

WORK  COMPLETED 

Since  the  start  of  this  current  award  the  following  tasks  have  been  completed: 

a.  Rewritten  the  hand-coded  TLM,  RPM,  and  ADM  algorithms  of  ROMS/TOMS  from  F77  to 
F90/F95  to  facilitate  multiple  levels  of  nesting,  parallelization  and  improve  computational 
efficiency.  The  TLM,  RPM  and  ADM  are  approximately  2.5,  2.6  and  2.8  times  more 
computationally  expensive  than  the  Nonlinear  Model  (NLM).  These  algorithms  were  updated  to 
the  latest  version  of  ROMS/TOMS  framework,  tested  and  distributed  to  selected  beta-testers 
around  the  world  on  May  15,  2006. 

b.  Parallelized  TLM,  RPM,  and  ADM  models  and  their  associated  drivers.  The  TLM  and  RPM 
can  be  run  in  either  shared-memory  (OpenMP)  or  distributed-memory  (MPI).  Currently,  the  ADM 
can  be  only  run  in  distributed-memory  because  its  hand-written  construction  violates  shared- 
memory  mutual  exclusion  rules  between  tiles.  The  TLM  and  RPM  have  a  parallel  structure 
identical  to  the  NLM.  The  ADM  required  adjoint  communication  exchanges  between  parallel 
domain  decomposition  tiles. 

c.  Developed  and  tested  the  strong  constraint,  incremental  4DVar  (IS4DVAR)  following  the 
approach  of  Courtier  et  al.  (1994).  The  modeling  of  the  background  error  covariance  uses  the 
generalized  diffusion  method  proposed  by  Weaver  and  Courtier  (2001).  The  free  parameters 
controlling  the  shape  of  the  Gaussian  correlation  for  each  state  variable  are  the  horizontal  and 
vertical  decorrelation  length-scales  and  the  diffusion  coefficients.  The  diffusion  operators  for  each 
state  variable  are  solved  explicitly  and  implicitly.  Since  the  oceanic  vertical  decorrelation  scales 
are  much  smaller  than  the  horizontal,  the  implicit  algorithm  is  preferable  and  cheaper.  It  is 
unconditionally  stable  for  any  vertical  convolution  time-step.  The  spatially  dependent 
normalization  coefficients  used  to  convert  the  covariance  matrix  into  a  correlation  matrix  are 
computed  using  the  exact  (expensive)  or  randomization  (cheaper)  methods  (Weaver  and  Courtier, 
2001).  These  normalization  coefficients  ensure  the  diagonal  elements  of  the  background  error 
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covariance  to  be  equal  to  unity.  After  extensive  testing  with  idealized  twin  experiments,  we  are 
now  working  with  realistic  data  assimilation  experiments  and  state  estimation  in  our  US  east  and 
west  coast,  Gulf  of  Mexico,  Caribbean  Sea,  and  East  Australia  Current  applications.  We  are 
assimilating  various  types  of  data  including  altimetry,  SST  from  satellites,  CTD,  XBT,  ADCP,  and 
gliders. 

d.  Developed  the  weak  constraint  4DVar  (W4DVAR)  single  and  multiple  drivers  using  the 
indirect  representer  approach  described  in  Chua  and  Bennett  (2001).  The  multiple  driver  option  is 
used  to  interface  with  the  IOM  framework  and  developed  under  separate  NSF/ITR  funding.  As  in 
IS4DVAR,  the  model  error  covariance  is  modeled  with  the  generalized  diffusion  operator. 
Currently,  we  are  using  the  W4DVAR  drivers  in  our  US  west  coast  and  Intra-Americas  Sea 
applications. 

e.  Developed  a  weak  constraint  4D-PSAS  (W4DPSAS)  driver.  The  PSAS  algorithm  is  similar  to 
W4DVAR  but  with  the  RPM  replaced  with  the  NLM.  That  is,  the  representer  functions  are  not 
explicitly  computed.  The  PSAS  acronym  is  misleading  but  it  is  retained  for  historical  reasons.  As 
in  W4DVAR,  the  minimization  is  in  observation  space.  Courtier  (1997)  shows  the  duality  between 
4DVar  and  PSAS.  Both  algorithms  produce  identical  results  if  the  measurement  functional  is  linear 
and  the  background  and  observation  error  covariance  are  the  same. 

RESULTS 

A  3D  baroclinic  coastal  upwelling  test  case  reveals  potentials  for  high  forecast  skill 

The  representer-based  weak  and  strong  constraint  4DVAR  ROMS  where  used  in  a  realistic  3D 
baroclinic  upwelling  system  with  complex  bottom  topography  (Di  Lorenzo  et  al.,  2006).  In  this 
example,  the  flow  field  is  nonlinear  and  characterized  by  mesoscale  activity  as  evident  from  the 
filamentary  and  cyclonic  structures  in  the  circulation  (Fig.  1).  Synthetic  observations  of  upper  ocean 
(0-45 0m)  temperatures  and  currents  with  a  high  resolution  sampling  array  (referred  to  as  the  HIRES 
sampling)  were  assimilated  over  a  10  day  window.  Both  the  strong  and  weak  constraint  inverse 
solutions  were  able  to  greatly  reduce  the  initial  error  variance  by  97%  and  80%  respectively  (Fig.  2). 
We  also  found  that  both  solutions  exhibit  relatively  high  forecast  skill  when  used  to  initialize  the 
nonlinear  model  at  the  end  of  the  assimilation  window.  Significant  forecast  skill  was  found  up  to  10-20 
days  after  the  last  observation  is  assimilated  and  is  higher  than  the  persistence  timescale  of  the  flow, 
which  for  this  upwelling  regime  is  less  than  5  days  (Fig.  2). 

The  same  experiment  was  repeated  using  an  observing  array  that  is  both  spatially  and  temporally 
aliased  (referred  to  as  the  COARSE  sampling).  For  this  case,  both  the  strong  and  weak  constraint 
posses  similar  levels  of  hindcast  and  forecast  skill,  although  the  weak  case  was  slightly  better.  During 
the  forecast,  the  skill  was  not  as  high  as  in  the  HIRES  case,  however  the  spatial  pattern  correlation  with 
the  true  is  still  very  high  when  compared  to  the  first  guess  (Fig  3).  These  results  suggest  that  the 
assimilation  platform  based  on  the  indirect  representer  method  with  inverse  ROMS  is  able  to  extract 
the  dynamically  active  information  from  the  observations  during  the  hindcast  window  and  generate  a 
good  initialization  for  the  forecast. 

The  application  of  the  assimilation  platform  to  this  highly  nonlinear  example  brought  to  our  attention 
one  important  aspect  of  the  methodology  that  required  further  consideration.  Specifically  we  found  that 
the  use  of  the  linearized  ROMS  model  (RP-ROMS)  in  an  iterative  approach  (the  outer  loop  in  the 
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representer  method)  does  not  necessarily  converge  to  the  solution  obtained  by  the  nonlinear  model. 

This  is  an  assumption  that  is  usually  made  in  the  representer  method.  Instead  we  found  that  RP-ROMS 
is  linearly  unstable  when  the  flow  field  is  very  nonlinear.  This  implies  that  in  the  assimilation  window 
the  corrections  provided  by  the  inverse  solution  need  also  to  damp  the  linear  instabilities  that  develop 
in  the  RP-ROMS.  In  our  current  work,  Di  Lorenzo  et  al.  (2006b)  different  formulations  of  the  outer 
loop  are  explored,  in  which  either  the  linear  instability  of  RP-ROMS  are  damped  by  including  a 
nudging  term  to  the  solution  of  the  previous  iteration  (Bennett,  2002), or  the  RP-ROMS  is  replaced 
with  NL-ROMS.  This  second  approach  is  similar  to  Courtier  (1997).  Preliminary  results  suggest  that 
these  approaches  greatly  improve  the  inverse  solutions. 

Real-time  forecast  for  the  Southern  California  Current  System  in  May  2006 

Applications  of  inverse  ROMS  using  real  ocean  observations  are  currently  being  performed  for  the 
Southern  California  Bight  (Di  Lorenzo  et  al.,  2006b)  and  for  the  Intra-Americas  Seas  (Powell  et  al., 
pers.  comm.).  Fig.  4  shows  our  first  real-time  forecast  for  the  oceanic  conditions  off  Southern 
California  during  May  2006.  The  forecast  was  initialized  using  the  ROMS  assimilation  platform  in 
combination  with  cruise  CTD  upper  ocean  data  for  the  month  of  April  2006.  The  forecast  was  posted 
on  the  web  (http://www.o3d.org/web/CalCOFI/april2006/)  on  April  30,  2005  and  was  used  to  guide  the 
biological  sampling  during  a  subsequent  cruise  in  May  2006.  Preliminary  analysis  of  the  forecast  skill 
reveals  that  the  model  was  able  to  predict  the  displacement  of  the  upwelling  front  (evident  in  Fig.  4) 
during  mid-may.  Further  analysis  are  currently  being  performed  to  fully  diagnose  the  model  prediction 
and  the  sources  of  dynamic  predictability  in  the  SCB. 

IMPACTS/APPLICATIONS 

The  newly  developed  ROMS/TOMS  adjoint-based  platforms  are  powerful  tools  for  ocean  prediction, 
adaptive  sampling,  and  understanding  the  underlying  circulation  dynamics. 

TRANSITIONS 

The  work  completed  here  will  be  part  of  the  ROMS/TOMS  utilities  that  will  be  freely  available  to  both 
research  and  operational  communities. 

RELATED  PROJECTS 

The  work  described  here  is  in  collaboration  with  Dr.  Andrew  Moore  at  University  California  at  Santa 
Cruz,  Drs.  Arthur  Miller  and  Bruce  Comuelle  at  the  Scripps  Institution  of  Oceanography,  and  Dr. 
Heman  Arango  at  Rutgers  University.  These  investigators  are  supported  by  the  following  grants: 

“ROMS/TOMS  Tangent  Linear  and  Adjoint  Models:  Testing  and  Applications”,  PI  Andrew  Moore, 
grant  number  N00014-01-1-0209. 

“ROMS/TOMS  Tangent  Linear  and  Adjoint  Models:  Testing  and  Applications”,  PI  Heman  Arango, 
grant  number  N00014-05-1-0366. 

“ROMS/TOMS  Tangent  Linear  and  Adjoint  Models:  Data  Assimilation  Tools  and  Techniques”,  Pis 
Arthur  Miller  and  Brace  Comuelle,  grant  number  N00014-99-1-0045. 
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Also  relevant  is  a  collaborative  project  involving  the  PI,  Dr.  Andrew  Moore  at  University  California  at 
Santa  Cruz,  Dr.  Heman  Arango  at  Rugers  University,  and  Dr.  Ralph  Milliff  at  Colorado  Research 
Associates: 

“Ocean  State  Estimation  and  Prediction  of  the  Intra-Americas  Sea”,  Pis:  Andrew  Moore,  Ralph  Milliff 
and  Heman  Arango,  grant  numbers:  N0014-05-M-0277,  N00014-05-M-0081,  N00014-05-0275. 
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Figure  1:  Maps  of  upper  ocean  temperature  ( 0-100  m)  for  hindcast  period  from 
DAY=0  (initial  condition)  to  DAY=10  (end  of  assimilation  window).  The  first 
column  represent  the  “true”  state,  the  second  column  is  the  basic  state  integration 
initialized  from  climatology,  the  third  column  is  the  result  from  the  weak  constraint 
assimilation  experiment  (Exp_weakH)  and  the  fourth  column  is  the  result  from  the 
strong  constraint  assimilation  experiment  (Exp_strongH) 
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Figure  2:  Hindcast  and  forecast  skill  scores  for  upper  ocean  temperature  ( 0-100  m) 
and  alongshore  velocity.  Skill  is  defined  based  on  the  RMS  difference  from  the 
truth  (see  section  5.6  in  text).  A  perfect  skill  value  is  1.  The  red  line  is  the  skill  of 
the  weak  constraint  solution  and  the  green  line  of  the  strong  constraint.  The  blue 
line  corresponds  to  the  skill  of persistence.  The  field  used  to  compute  the 
persistence  skill  is  computed  by  taking  the  available  observations  and  performing 
an  objective  mapping  to  interpolate  at  the  location  where  observations  are  not 
available.  In  the  forecast  window  (DAY  10-30)  no  observations  are  used  to 
constrain  the  model  trajectory.  The  upper  row  shows  result  using  the  HIRES 
sampling  array  and  the  lower  row  the  COARSE  array. 
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Figure  3:  Maps  of  upper  ocean  temperature  (0-100  m)  in  the  forecast  window  from 
DAY=14  (4  days  after  the  end  of  the  assimilation  window)  to  DAY=26.  The  first  column 
represent  the  “true”  state,  the  second  column  is  the  basic  state  integration  initialized  with 
climatology,  the  third  column  is  the  result  from  the  weak  constraint  assimilation 
experiment  (Exp_weakC)  and  the  fourth  column  is  the  result  from  the  strong  constraint 
assimilation  experiment  ( Exp_strongC ').  The  forecast  are  performed  using  the  nonlinear 
model  initialized  with  the  solutions  of  the  assimilation  experiments  at  DAY=10. 
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Figure  4:  An  image  of  the  on-line  real  time  ocean  forecast  for  May  2006.  The 
forecast  was  initialized  using  4D  variaitonal  data  assimilation  with  the  inverse 
Regional  Ocean  Modleing  System  (ROMS).  During  the  assimilation  time  window, 
TS  data  from  the  April  2006  CalCOFI  cruise  where  assimilated  to  initialize  the 
model  field  for  the  forecast.  In  the  forecast,  the  model  was  forced  with  winds 
provided  by  the  ECPC  regional  atmospheric  model  forecast.  This  forecast  was  used 
to  guide  biological  oceanographers  during  a  cruise  in  May  2006. 
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